Dispersion properties of electrostatic oscillations in quantum plasmas 
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We present a derivation of the dispersion relation for electrostatic oscillations (ESOs) in a zero 
temperature quantum plasma. In the latter, degenerate electrons are governed by the Wigner equa- 
tion, while non-degenerate ions follow the classical fluid equations. The Poisson equation determines 
the electrostatic wave potential. We consider parameters ranging from semiconductor plasmas to 
metallic plasmas and electron densities of compressed matter such as in laser-compression schemes 
and dense astrophysical objects. Due to the wave diffraction caused by overlapping electron wave 
function due to the Heisenberg uncertainty principle in dense plasmas, we have possibility of Landau 
damping of the high-frequency electron plasma oscillations (EPOs) at large enough wavenumbers. 
The exact dispersion relations for the EPOs are solved numerically and compared to the ones ob- 
tained by using approximate formulas for the electron susceptibility in the high- and low-frequency 



I. INTRODUCTION 



The field of quantum p lasma ph y sics is becoming of increasing; cur r ent interest Bonitz et al. I (120031k Crouseilles 



20081 ): iManrredil (I2005D ; fSerbeto I (I2008D : IShaikh and Shukla ] (l2007h : IShukla and Eliasson I (|2006h : IShukla I (12009)]. 



motivated by its potential applications in modern technology (e.g. metallic and semiconductor nanostructures-such 
as metallic nanoparticles, metal clusters, thin metal films, spintronics, nanotubes, quantum well and quantum dots, 
nano-plasmonic devices, quantum x-ray free-electron lasers, etc.). In dense quantum plasmas and in the Fermi gas 
of metals, the number densities of degenerate electrons are extremely high so that their wave functions overlap, and 
they therefore electrons obey the Fe rmi-Dirac stat i stics. The collective osc il lations i n quantum plasmas have been 
studied by s e veral authors in the past [Bohm I (|19520 ; iBohm and Pines 1 (|l953l ): lFerrell (|l957D ; lKlimontovich and Silin I 
(jl952lll96ll ): lPmesl(|l96ll) ] with appl ications to the Fe r mi pla s mas in metals and semiconductors, and to electrostatic 
oscillations in quantum pair plasmas [Mendonca et all ( 20081 )]. IWatanabel ( 19561 ) studied experimentally the Bohm- 
Pines dispersion relation of the electron plasma oscillations by measuring the energy loss of electrons by the excitation 
of collective modes in metals. The Fermi degenerate dense plasma may also arise when a pellet o f hydrogen is 
compr e ssed to many t i mes the solid density i n the fast igniti o n scen a rio for inertia! confinement fusion Azechie£ al. I 
(119911) : lAzechi et al7\ (|2006l) : iLindl I (| 19951) : ISon and Fischl (|2005l) : iTabak et al. 1 (11994 I2005D ]. Since there is an 
impressive developments in the field of short pulse petawatt laser technology, it is highly likely that such plasma 
conditions can be achieved by intense laser pulse compression using powerful x-ray pulses. Here ultrafast x-ray 
Thomson scattering techniques can be used to measure the features of laser enhanced plasma lines, which will, in 
turn, give invaluable informations regarding the equation of state of s h ock compressed d ense matters. Recently, 
spectrally resolved x-ray scattering measurements [Kritcher et al. I ( 20081 ): iLee et al~\ (|2009[ )] have been performed in 
dense plasmas allowing accurate measurements of the elect ron velocity dis t ributi on function, temperature, ionization 
state, and of plasmons in the warm dense matter regime [Glenzer et al. I (2007)]. This novel technique promises to 
access the degenerate, the closely coupled, and the ideal plasma regime, making it possible to investigate extremely 
dense states of matter, such as the inertial confinement fusion fuel during compression, reaching super-solid densities. 

In this paper, we present a study of the dispersion properties of electrostatic oscillations in a dense quantum 
plasma, by employing the Wigner-Poisson model. We point out the differences between different regimes comprising 
the relatively low density regime of semiconductor plasmas, and the higher density regimes corresponding to metallic 
electron densities and laser compressed plasmas, as well as plasmas in dense astrophysical objects such as white dwarf 
stars. 



II. DERIVATION OF THE DISPERSION RELATION FOR THE WIGNER-POISSON SYSTEM 



We here present a derivation of the dispersion relation for electrostatic waves in a degenerate quantum plasma. The 
electron dynamics is governed by the Wigner equation 
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where the electrostatic potential is given by the Poisson equation 
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Here e is the magnitude of the electron charge, m e is the electron mass, h is the Planck constant divided by 2tt, and 
6o is the permittivity of free space. Furthermore, /o and no denote the equilibrium electron distribution function and 
the electron number density, respectively, while /i, 4>\ and rin denote the perturbed electron distribution function, 
the electrostatic potential, and the ion number density, respectively. 

Assuming that f\, 4>\ and tin are proportional to exp(— itot + ik • x), where u> is the frequency and k is the wave 
vector, we obtain from Eq. (fTJ) and ©, respectively, 
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Since ions are non-degenerate in quantum plasmas, we have for u> 3> kVn 
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is the ion susceptibility, Vy$ is the ion thermal speed, and uj pi is the ion plasma frequency. 
Rewriting (J3J) as 
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and performing the integration over A space, we have 



d 3 Xd 3 v' 



i— (v — v') ■ A — zk • A/ 



/2]}/ o (v')0i(w,k), 



(7) 



(w - v k)A = — jS 



m e , , k 
X (v - v) -2 



/o(v')<M^k), 



X (V ^ v)+ 2 



(8) 



where 5 is the Dirac delta function. Now, the integration can be performed over v' space, obtaining the result 
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Eliminating na and /i in (J4j> with the help of (J5J and ((9]), we obtain the dispersion relation 

1 + Xe + Xi = 0, 

where the ion susceptibility is given by ((6|) and the electron susceptibility is given by 
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Suitable changes of variables in the two terms in square brackets in Eq. (jlip now give 
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This expression was also been derived by iBohm and Pines ] (119531) by using a series of canonical t ransformation s of 
the Hamiltonian of the system [see for example the dispersion relation (57) in their paper], and bv lFerrell (|l957h by 
using the method of self-consistent fields. 

We now choose a coordinate system such that the x axis is aligned with the wave vector k. Then, (fT3]) takes the 
form 
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We next consider a dense plasma with degenerate electrons in the zero temperature limit. Then, the background 
distribution function takes the simple form 
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where Vf £ = i^Fe/me) 1 ^ 2 is the speed of an electron on the Fermi surface, and £fc = (3ir 2 no) 2 / 3 h 2 / (2m e ) is the 
Fermi energy. The integration in (|14[) can be performed over velocity space perpendicular to u x , using cylindrical 
coordinate in u y and u z , obtaining the result 
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It is interesting to note that the distribution, which is flat-topped in three-dimensions becomes parabola-shaped in 
the remaining velocity dimension after the integration over the two perpendicular velocity dimensions. Hence, the 
electron distribution function Fq(u x ) in (fl7| may support Landau damping if the pole of the denominator in (fl6|) falls 
into the range of negative slope of Fq(u x ) in velocity space. Equation (|16[) can be written as 
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Performing the integration over velocity space, we have from (fl8 
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In the derivation of (|19[) . we have assumed that the waves are only weakly damped, so that, when integrating over 
poles, only the principal parts of the integrals are kept. In the limit hk/m e — > 0, we have from (|19[) 
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where it holds that lo is real and uj/k > Vfc- Hence in this "semi-classical" limit, we do not have Landau damping. 
[We call ([2H)l "semi-classical", since it can be derived from the Vlasov equation for electrons using the flat-topped 
background electron distribution function given by (|15p .] 



III. ELECTRON OSCILLATIONS 



We here consider high-frequency (uj 3> Lo p i) oscillations so that ^ <C 1 in ([BJ. Then, expanding (jT9|) for small 
wavenumbers up to terms containing fc 4 , we obtain from (|10| the dispersion relation 
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where a = (48/ 17 '5)m 2 e V^ e /h 2 ujl e w 2.000(n o ao) 1/3 and a Q = h 2 /m e e 2 w 53 x 10" 10 cm is the Bohr radius. For 
a typical metal such as gold, which has a free electron number density of no = 5.9 x 10 22 cm -3 , we would have 
a w 0.4. For the free electron density in semiconductors, which is many orders of magnitude less than in metals, 
a is much smaller and can safely be dropped compared to unity. However, for electron plasma oscillations in dense 
m atters, a could be l arger than unity. It should be noted that the term proportional to a in (1211) w as not discussed 
by iBohm and Pines I |l953) and others, but was, however, obtained and discussed by iFerreil ( 19571 ) in his study of 
collective electron oscillations in metals [see Eq. (10) in his paper, where, in his notation, it should be (Av 2 /v 2 ) 2 = 
12/175]. 

It is interesting to note that (|19p may admit Landau damping above a certain critical wavenumber, k > k cr and 
corresponding frequency w > u> cr . This occurs if the denominator in the integral of (|18p vanishes within the integration 
limits u x = ±V Fe . For the critical wavenumber and frequency, we have 
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Inserting this expression into (|19p we note that the term involving the logarithm on the second line of (|19p vanishes, 
and we obtain the critical wavenumber k = k cr from 
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A careful examination of the dispersion relation for fc > k cr , should involve Landau contours to correctly take into 
account Landau damping. Here we are interested in low- and high-frequency waves in the weakly damped regime and 
have postponed the investigation of Landau damping of the system to future studies. 

IV. ION OSCILLATIONS 

In a quantum plasma system composed of mobile ions and inertialess electrons, we have possibility of low-phase speed 
(in comparison with the Fermi electron thermal speed) ion-acoustic- like oscillations. For low- frequency (u> <C kVp e ) 
waves, we have from (|19[) 
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For small wavenumbers hk <C m e VF e , we have the approximate electron susceptibility, up to terms containing factors 
of fc 4 , 
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we employ (j^J) to obtain the frequency of ion acoustic waves as 



(25) 



(26) 



to 



-"pi 



(1 + Xe) 1 / 2 ' 

with \e given by (|24[) or (|25p . For Xe given by the approximate expression (|25p . we have 
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where C s — \J m e Vp e /3m.i is the Fermi ion acoustic speed. We note that 



lu — > ujpi as k — » co. 



V. DENSITY REGIMES OF THE SYSTEM 

We note that there is a critical density parameter in the system. When the inter-particle distance is smaller than 
the Bohr radius, then the quantum statistical pressure dominates the wave dynamics, while in the opposite case, the 
quantum tunneling effects become important when the wavelength is comparable to the inter-particle distance. This 
can be seen by normalize the system such that u>/oj pe = Q and fcVp e /w pe = K. Then, Eq. (|19[) takes the form 
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In the scaled variables, (l20l) and (12TT) take the form 
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respectively. We note that the limit hk/m e — > to obtain (f20f from (|19p corresponds to /3 — > oo to obtain (|3~lj) 
from (|29| . For j3 > 1, the quantum statistical pressure dominates, while for (3 <C 1, the quantum tunneling effects 
dominate. Considering the value of j3 for different physical systems, we note that (3 — 0.1 corresponds to relatively 
low density degenerate plasma such as in semiconductors, while (3 = 1 corresponds to typical free electron densities in 
metals. The high density case /? = 10 corresponds to high density matter which may be obtained in laser compression 
schemes or which exist in white dwarf stars. Even though we have formally considered the semiclassical limit (3 — > oo, 
it should be kept in mind that an upper limit for the validity of our theory is when the electron density becomes high 
enough that the Fermi speed Vp e becomes comparable with the speed of light. In this limit, the inter-particle distance 
n" 1 / 3 approaches the Compton length Ac = 2nh/m e c w 2.4 x 10~ 12 , and we have an electron number density of the 
order 10 35 m~ 3 , corresponding to (3 w 27. For larger values of (3, the equilibrium equation of state for the electrons 



ChandrasekhaTl dl935h ] change from P = (2/5)£ F no(n e /n ) 5/3 to P = (3/7r) 1 / 3 (47rc/i/8) n e ■ For this case, we need 



to include relativistic effects in the electron susceptibility. 

In the normalized variables, the condition (|23p for the critical wavenumber for the limit between undamped and 
Landau damped high-frequency waves is given by 
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and the normalized critical frequency is obtained from 
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where K cr — k cr VFe/^pe and fl cr — uj cr /uj pe . 

Finally, for the low-frequency case, the electron susceptibilities (|24p is in the normalized variables given by 



2K 2 



and the expression for small wavenumbcrs (|25[) is given by 



1 - 


K 

2/3 




1 + 


K 

2/3 





K 2 + TiK 4 / '12 (3 2 



(36) 



VI. NUMERICAL RESULTS 



In Figs. 1 and 2, we show dispersion curves for the high-frequency (lu 3> uj p i) waves for different values of /3, 
obtained from the solutions of the dispersion relation (fTT)]) . by using the electron susceptibility (HHJ), as well as the 
expansion (|2"Tj) and the limiting semi-classical case ([2ll|) . We have also indicated the border between undamped and 
Landau damped waves, obtained from (f2"2")) and (|2"3"|) . We note that the dispersion curve for the semiclassical case 
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FIG. 1: Dispersion curves (w versus k) for different values of f3 = m e Vp e /hujp e . The solid curves show solutions of (|10[) using 
the exact susceptibility (|19p , the dashed curves show the expanded solution (|2ip . and the dash-dotted curve shows solutions 
of (|10[) using the "semi-classical" electron susceptibility (|20[) . The dotted curve indicates the border between undamped waves 
and Landau damped waves, given by (3.2) and (3.3) [or (5.5) and (5.6)]. 
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FIG. 2: Closeup of the dispersion curves in Fig. 1. 



in Fig 1 always lies in the undamped regime, below the border between undamped and Landau damped waves. For 
the undamped waves, the expansion (|2ip approximates the exact dispersion relation within a few percent, and can 
therefore be used instead of p^|) for most cases. This holds especially for small wavenumbers, as can be seen in the 
closeup in Fig. 2. 

The dispersion curves for the low-frequency ion-acoustic oscillations are plotted in Fig. 3, were we have depicted 
the u> in f2"7) for \e given by the exact expression (pM| and the approximate expansion ([23]) . and for different values 
of (3. We note that the dispersion curves show agreement at small wavenumbers but deviate significantly for larger 
wavenumbers in the case /3 = 0.1, while the agreement is better for (3 — 1 and excellent for (3 = 10. It should be kept 
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FIG. 3: Dispersion curves (lu versus k) for the low-frequency ion oscillations, for different values of (3 — m e Vp e /hujp e . The solid 
curves show the wave frequency given by Q270 using the low- frequency electron susceptibility (|24p , while the dashed curves uses 
the approximate electron susceptibility (|25[1 . For /3 = 10 the solid and dashed curves are indistinguishable. 

in mind that when fik ~ m e VF e then the wavelength of the oscillations is comparable to the inter-particle distance, 
and there will be corrections due to the discrete nature of the ion background. Thus, the theory is not valid for 
oscillations with Hk 3> m e VF e . 

VII. CONCLUSIONS 

In this paper, we have here studied the dispersion properties of electrostatic oscillations in quantum plasmas 
for different parameters ranging from semiconductor plasmas, to typical metallic electron densities and to densities 
corresponding to compressed matter and dense astrophysical objects. We have derived a simplified expansion that 
accurately approximates the exact dispersion relation for small wavenumbers. The possibility of Landau damping 
due to quantum tunneling effects at large wavenumbers has also been discussed, and conditions for Landau damping 
has been derived. The present results should be useful in understanding the salient features of electrostatic plasma 
oscillations in dense plasmas with degenerate electrons. The latter are encountered in metals, in highly compressed 
intense laser-solid density plasma experiments, and in compact astrophysical objects (e.g. interior of white dwarf 
stars). 
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